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Abstract 

Review of recent works devoted to the temporal and spatial variation of the fundamental con- 
stants and dependence of the fundamental constants on the gravitational potential (violation of 
local position invariance) is presented. We discuss the variation of the fine structure constant 
a = e 2 /he, strong interaction and fundamental masses (Higgs vacuum), e.g. the electron-to-proton 
mass ratio \x = m e /M p or X e = m e / Aqcd and X q = m q / 'Aqcd- We also present new results from 
Big Bang nucleosynthesis and Oklo natural nuclear reactor data and propose new measurements 
of enhanced effects in atoms, nuclei and molecules, both in quasar and laboratory spectra. 

PACS numbers: 



I. INTRODUCTION 



Theories unifying gravity with other interactions suggest temporal and spatial variation 
of the fundamental "constants" in expanding Universe (see e.g. review The spatial 
variation can explain fine tuning of the fundamental constants which allows humans (and 
any life) to appear. We appeared in the area of the Universe where the values of the 
fundamental constants are consistent with our existence. The fundamental constants may 
be slightly different near massive bodies (see e.g. review (2|). There are some hints for the 
variation of different fundamental constants in quasar absorption spectra 
and Big Bang nucleosynthesis jsl, Q data. However, a majority of publications report limits 
on the variations (see e.g. recent reviews [ill . [T" 



The hypothetical unification of all interactions imp 



ies that variations of different fun- 



damental constants may be related Q, Q, Q, 12 1- We can only detect variation of 
dimensionless fundamental constants. We will discuss variation of the fine structure con- 
stant a and dimensionless ratios X e = m e / Aqcd and X q = m q / Aqcd where m e and m q 
are the electron and quark masses, and Aqcd is the quantum chromodynamics (QCD) scale 
(defined as the position of the Landau pole in the logarithm for the running strong coupling 
constant, a s {r) ~ 1/ln (Aqcdt /he)). The proton mass m p is proportional to Aqcd, there- 
fore the relative variation of fi = m e /M p is equal to the relative variation of X e = m e / Aqcd 
(if we neglect a small contribution of quark masses (m q ~ 5 MeV) to the proton mass, 
m p = 938 MeV). In the Standard model electron and quark masses are proportional to the 
vacuum expectation value of the Higgs field. 

A simple estimate of the relations between the variations of different fundamental con- 
stants may be obtained using the idea of Grand Unification. The strong (i=3), and elec- 
troweak (i=l,2) inverse coupling constants have the following dependence on the scale v and 
normalization point u : 

arV) = a* 7 Vo) + biln(u/v ) (1) 

In the Standard Model, 2nbi = 41/10, —19/6, —7; the electromagnetic or 1 = (h/tyai 1 + a^ 1 
and the strong a s = 03. In the Grand Unification Theories (GUT) all coupling constants are 
equal at the unification scale, a^o) = &gut- If we assume that oigut varies, then Eq. ([1]) 
gives us the same shifts for all inverse couplings: 

Sa^ 1 = 5a 2 1 = Sa^ 1 = ^a^c/T • (2) 



We see that the variation of the strong interaction constant 0:3(1/) at low energy v 
is much larger than the variation of the elecromagnetic constant a, since Sa^/a^ = 

(a 3 /a;i ) 2)5ai ) 2/ai,2 and «3 > «i,2- 

The variation of mj Aqcd can be estimated from the definition of Aqcd- The running of 
a s near the electroweak scale is given by 



a s {v) 1 ps b s ln(i> / Aqcd) 

Let us take z/ = m z where m z is the Z-boson mass. The variation of eq. 
above give 



(3) 



and relations 



S(m z /A QC D) 



5a s {m z ) 



C 5mm, 



(4) 



(m z /A QCD ) b s a s (m z ) a s {m z ) a{m z ) a{m z ) 
The value of the constant C here depends on the model used. However, the enhancement 
1/a ~ 100 should make the factor C/a large. Note that m z (as well as m e and m q ) is 
proportional to the Higgs vacuum expectation value. 

If this estimate is correct, the variation in X e>q = vn e ^j 'Aqcd or \i — m e /M p may be 
easier to detect than the variation in a. The cosmological variation of m q j Aqcd can be 
extracted from the big bang nucleosynthesis (BBN), quasar absorption spectra and Oklo 
natural nuclear reactor data 19, 2o| . For example, the factor of three disagreement 



between the calculations and measurements of the BBN abundance of 7 Li may, in principle, 

on 

be explained by the variation of m q / Aqcd at the level of ~ 10 3 — 10 2 |9j, llOj (see also recent 
work 2 IIP. The claim of the variation of the fundamental constants based on the Oklo data 



in Ref . 



22J is not confirmed by recent studies 23|, |2J, [25( which give a stringent limit on the 
possible variation of the resonance in 150 Sm during the last two billion years. The search 



for the variation of m e / Aq C d using t 



Ref. [8|] but zero results in Refs. [26 



re quasar absorption spectra gave a non-zero result in 
present time variation of m P n /A n cD can 



27 



The 



be extracted from comparison of different atomic 
New enhanced effects have been proposed. 



291 ] , molecular or nuclear [3l|, |32j clocks. 



II. BIG BANG NUCLEOSYNTHESIS 



The result of our work Ref. [9j suggested that a reduced deuteron binding energy of 
AQ/Q = -0.019 ± 0.005 would yield a better fit to observational data (the WMAP value 



3 



of barion-to-photon ratio?? and measured 2 H, 4 He, and 7 Li abundances) for Big Bang Nu 



20j we obtained in Ref. {9] an estimate of the strange 



cleosynthesis. Using our calculations 
quark mass variation. 

Recently Dent, Stern, and Wetterich 2l| calculated the sensitivity of BBN abundances 
for 2 H, 4 He and 7 Li to the variation of binding energies of 2 ' 3 H, 3,4 He, 6,7 Li and 7 Be in a 



linear approximation. In the works 



10 



33 



3, 



35] we calculated dependence of these binding 



energies on the light quark mass variation and estimated the sensitivity of BBN yields to 
variation of the quark mass. Then we used the observational data to obtain the following 
equations for 2 H, 4 He and 7 Li lol |: 

2.8 ±0.4 



1 + 7.7x 
1 _ o.95x = 
1 -50a; 



2.61 ±0.04 
0.249 ± 0.009 
0.2478 ±0.0002 
1.5 ±0.5 



1.07 ±0.15 , 

= 1.005 ±0.036 , 
0.33 ±0.11 , 



(5) 

(6) 
(7) 

4.5 ±0.4 "' ' v ; 

where x = 5X q /X q . These equations yield 3 consistent values of x: 0.009 ± 0.019, —0.005 ± 

0.038 and 0.013 ± 0.002. The statistically weighted average of 5X q /X q = 0.013 ± 0.002 

is dominated by the 7 Li data. A more accurate calculation should take into account the 

effect of the 8 Be binding energy variation (which is not calculated in Ref. 2l|), the variation 

of the virtual 1 So(?^p) level, and non-linear corrections in x which are important for 7 Li. 

Allowing for the theoretical uncertanties we should understand this BBN result as 5X q /X q = 

K ■ (0.013 ± 0.002) where K ~ 1, where the expected accuracy in K is about a factor of 2. 

Note that here we neglected effects of the strange quark mass variation. A rough estimate of 

these effects on BBN due to the deuteron binding energy variation was made in Refs. 0, Q] . 



III. OKLO NATURAL NUCLEAR REACTOR 



The results from Oklo natural nuclear reactor are based on the measurement of the 
position of very low energy resonance (E r = 0.1 eV) in neutron capture by 149 Sm nucleus. 
The estimate of the shift of this resonance induced by the variation of a have been done 
long time ago in works 



variation of m q j 'Aq C d 



36, 

at 



37|. Recently we performed a rough estimate of the effect of the 
2o|. The final result is 



5E r « 10 6 eU(^ - 10^ 
a X n 



+ 100- 



.SXs 



(8) 



where X q = 
mass. Refs. 



m, 



23 



/Aopii X s = m s /A Q cD, 



24 



(m„ + md)/2 and m s is the strange quark 



251 ] found that \8E r \ < 0.1 eV. This gives us a limit 



<5a 
0.01 — 

a 



0.1 9 



+ 



X q X s 



< io- 



cs) 



The contribution of the a variation in this equation is very small and should be neglected 
since the accuracy of the calculation of the main term is low. Thus, the Oklo data can not 
give any limit on the variation of a. Assuming linear time dependence during last 2 billion 
years we obtain an estimate \X S /X S \ < 10~ 18 yr -1 . 

IV. OPTICAL ATOMIC SPECTRA 

A. Comparison of quasar absorption spectra with laboratory spectra 



To perform measurements of a variation by comparison of cosmic and laboratory optical 
spectra we developed a new approach 



38 



39] which improves the sensitivity to a variation 



of a by more than an order of magnitude. The relative value of any relativistic corrections 
to atomic transition frequencies is proportional to a 2 . These corrections can exceed the 
fine structure interval between the excited levels by an order of magnitude (for example, 
an s-wave electron does not have the spin-orbit splitting but it has the maximal relativistic 
correction to energy). The relativistic corrections vary very strongly from atom to atom 
and can have opposite signs in different transitions (for example, in s-p and d-p transitions). 
Thus, any variation of a could be revealed by comparing different transitions in different 
atoms in cosmic and laboratory spectra. 

This method provides an order of magnitude precision gain compared to measurements 
of the fine structure interval. Relativistic many-body calculations are used to reveal the 



dependence of atomic 
absorption spectra 38 



requencies 



39 



42 



on a for a range of atomic species observed in quasar 



43 



44 , |45| (a 2004 summary may be found in Ref. [46| ) . 
It is convenient to present results for the transition frequencies as functions of a 2 in the form 



U! 



uj + qx, 



(10) 



where x = (^) 2 — 1 ~ 222 anc l Uq j s a laboratory frequency of a particular transition. 



We stress that the second term contributes only if a deviates from the laboratory value 
a . We performed accurate many-body calculations of the coefficients q for all transtions of 



astrophysical interest (strong El transtions from the ground state) in Mg, Mg II, Fe II, Fe 
I, Cr II, Ni II, Al II, Al III, Si II, Zn II, Mn II (and many other atoms and ions which have 
not been used in the quasar measurements yet because of the absence of accurate laboratory 
wavelenths - see 46]). It is very important that this set of transtions contains three large 
classes : positive shifters (large positive coefficients q > 1000 cm -1 ), negative shifters (large 
negative coefficients q < —1000 cm -1 ) and anchor lines with small values of q. This gives 
us an excellent control of systematic errors since systematic effects do not "know" about 
sign and magnitude of q. Comparison of cosmic frequencies uj and laboratory frequencies u 
allows us to measure — . 

a 

Three independent samples of data contaning 143 absorption systems spread over red 
shift range 0.2 < z < 4.2. The fit of the data gives |5( is ^ = (-0.543 ± 0.116) x 10~ 5 . 
If one assumes the linear dependence of a on time, the fit of the data gives dlna/dt = 
(6.40 ± 1.35) x 10~ 16 per year (over time interval about 12 billion years). A very extensive 
search for possible systematic errors has shown that known systematic effects can not explain 
the result (It is still not completely excluded that the effect may be imitated by a large 
change of abundances of isotopes during last 10 billion years. We have checked that different 
isotopic abundances for any single element can not imitate the observed effect. It may be 
an improbable "conspiracy" of several elements) . 



Recently our method and calculations 38l . 



groups 



39 



42 



43 



44j were used by two other 



471. 1481 . 491. However, they have not detected any variation of a. Recently the 

tt n q n 

results of [47J were questioned in Refs. |6j, ITJ • Re-analysis of Ref . [47J data revealed flawed 
parameter estimation methods. The authors of js, 2] claim that the same spectral data fitted 
more accurately give % = (-0.64 ± 0.36) x 10~ 5 (instead of — = (-0.06 ± 0.06) x 10~ 5 in 



Ref. 14 



However, even thisjreyised result may require further revision. 

Xeck te^ 



47 



48 



escope which 



491 ] are based 



Note that the results of are based on the data from the 

is located in the Northen hemisphere (Hawaii). The results of 
on the data from the different telescope (VLT) located in the Southern hemisphere (Chile). 
Therefore, some difference in the results may appear due to the spatial variation of a. 

Using opportunity I would like to ask for new, more accurate laboratory measurements of 
UV transition frequencies which have been observed in the quasar absorption spectra. The 
"shopping list" is presented in [46|. We also need the laboratory measurements of isotopic 
shifts - see (46[. We have performed very complicated calculations of these isotopic shifts 



6 



50, |5l|, |52|, |53j, l54(. However, the accuracy of these calculations in atoms and ions with 
open d-shell (like Fe II, Ni II, Cr II, Mn II, Ti II) may be very low. The measurements for 
at list few lines are needed to test these calculations. These measurements would be very 
important for a study of evolution of isotope abundances in the Universe, to exclude the 
systematic effects in the search for a variation and to test models of nuclear reactions in 
stars and supernovi. 

A comparison of the hyperfine transition in atomic hydrogen with optical transitions 
in ions, was done in Refs. 26j, |27j. This method allows one to study time- variation of 



the parameter F = a 2 g p fi, where g p is proton ^-factor. Analysis of 9 quasar spectra with 
redshifts 0.23 < z < 2.35 gave 



SF/F= (6.3 ±9.9) x 10" 



F/F = (-6 ± 12) x 10~ 16 yr" 1 



(11) 
(12) 



B. Optical atomic clocks 

Optical clocks also include transitions which have positive, negative or small constri- 
butions of the relativistic corrections to frequencies. We used the same methods o 
relativistic many-body calculations to calculate the dependence on a J39 



A 2004 summary of the results for the coefficients q is presented in [58]. The coefficients q 
for optical clock transitions may be substantially larger than in cosmic transitions since the 
clock transitions are often in heavy atoms (Hg II, Yb II, Yb III, etc.) while cosmic spectra 
contain mostly light atoms lines (Z < 33). The relativistic effects are proporitional to Z 2 a 2 . 

V. ENHANCED EFFECTS OF a VARIATION IN ATOMS 



An enhancement of the relative effect of a variation can be obtained in transition between 
the almost degenerate levels in Dy atom 



39, 



40 



571 ] . These levels move in opposite directions 



if a varies. The relative variation may be presented as Su/u = KSa/a where the coefficient 
K exceeds 10 8 . Specific values of K = 2q/u are different for different hyperfine components 



and isotopes which have different u; q = 30, 000 cm 



-i 



10 4 cm 1 . An experiment is 



currently underway to place limits on a variation using this transition 



59, 



60j. The current 
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limit is a/ a = (—2.7 ± 2.6) x 10~ 15 yr" 1 . Unfortunately, one of the levels has quite a large 
linewidth and this limits the accuracy. 

Several enhanced effects of a variation in atoms have been calculated in 6l|, |62|. 



VI. ENHANCED EFFECT OF VARIATION OF a AND STRONG INTERACTION 
IN UV TRANSITION OF 22 9TH NUCLEUS (NUCLEAR CLOCK) 



A very narrow level (7.6 ± 0.5) eV above the ground state exists in 229 Th nucleus 63]. 
The position of this level was determined from the energy differences of many high-energy 7- 
transitions to the ground and excited states. The subtraction produces the large uncertainty 
in the position of the 7.6 eV excited state. The width of this level is estimated to be about 
10 -4 Hz |64j. This would explain why it is so hard to find the direct radiation in this very 



65 1 



weak transition. However, the search for the direct radiation continues 

The 229 Th transition is very narrow and can be investigated with laser spectroscopy. This 
makes 229 Th a possible reference for an optical clock of very high accuracy, and opens a new 



possibility for a laboratory search for the varitation of the fundamental constants 



As it is shown in Ref. 



32 1 



3JJ there is an additional very important advantage. According 
to Ref. |3j| the relative effects of variation of a and m q / Aqcd are enhanced by 5 orders of 
magnitude. This estimate has been confirmed by the recent calculation [^j and preliminary 
results of our new calculations. The accurate results of the calculations will be published 
soon. A rough estimate for the relative variation of the 229 Th transition frequency is 

Sou K . 5a SX n . 

— « 10 5 0.1- + — 1) (13) 
u a X q 

where X q = m q / Aqcd- Therefore, the Th experiment would have the potential of improv- 
ing the sensitivity to temporal variation of the fundamental constants by many orders of 
magnitude. Indeed, we obtain the following energy shift in 7.6 eV 229 Th transition: 

5uj « — q -MeV (14) 

This corresponds to the frequency shift 8u ~ 3 • 10 20 5X q / X q Hz. The width of this transition 
is 10~ 4 Hz so one may hope to get the sensitivity to the variation of X q about 10 -24 per 
year. This is 10 10 times better than the current atomic clock limit on the variation of X q , 
~ 10 -14 per year. 



S 



Note that there are other narrow low-energy levels in nuclei, e.g. 76 eV level in 235 JJ 
with the 26.6 minutes lifetime (see e.g. 32J). One may expect a similar enhancement there. 
Unfortunetely, this level can not be reached with usual lasers. In principle, it may be 
investigated using a free-electron laser or synchrotron radiation. However, the accuracy of 
the frequency measurements is much lower in this case. 



VII. ATOMIC MICROWAVE CLOCKS 



Hyperfine microwave transitions may be used to search for a-variation 



671 ] . Karshenboim 



681 ] has pointed out that measurements of ratios of hyperfine structure intervals in different 



atoms are also sensitive to variations in nuclear magnetic moments. However, the magnetic 
moments are not the fundamental parameters and can not be directly compared with any 
theory of the variations. Atomic and nuclear calculations are needed for the interpretation 



of t 
40, 



re measurements. We have performed both atomic calculations of a dependence 39, 



41 



n 



56 



57 



and nuclear calculations of X n 



m, 



,/Aqcd dependence [29| (see 



also [35j) for all microwave transitions of current experimental interest including hyperfine 
transitions in 133 Cs, 87 Rb, m Yb+, 199 Hg+, m Cd, 129 Xe, 139 La, 1 H, 2 H and 3 He. The results 
for the dependence of the transition frequencies on variation of a, X e = m e /A.QCD and 
X q = m q / Aqcd are presented in Ref.[29j (see the final results in the Table IV of Ref . |29| ) . 
Also, one can find there experimental limits on these variations which follow from the recent 
measurements. The accuracy is approaching 10~ 15 per year. This may be compared to the 
sensitivity ~ 1CT 5 — 10~ 6 per 10 10 years obtained using the quasar absorption spectra. 

According to Ref. 29| the frequency ratio Y of the 282-nm 199 Hg + optical clock transi- 
tion to the ground state hyperfine transition in 133 Cs has the following dependence on the 
fundamental constants: 

Y/Y = -6a/a - fi/fi - 0MX q /X q (15) 

In the work 30] this ratio has been measured: Y/Y = (0.37 ±0.39) x 10~ 15 yr . Assuming 
linear time dependence we obtained the quasar result [28| fi/fi = XJX e = (1 ± 3) x 
10~ 16 yr -1 . A combination of this result and the atomic clock result 30| for Y gives the 



best limt on the variation of a: a/a = (—0.8 ± 0.8) x 10 16 yr 1 . Here we neglected the 
small (~ 1%) contribution of X q . 
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VIII. ENHANCEMENT OF VARIATION OF FUNDAMENTAL CONSTANTS IN 
ULTRACOLD ATOM AND MOLECULE SYSTEMS NEAR FESHBACH RESO- 
NANCES 



Scattering length A, which can be measured in Bose-Einstein condensate and Feshbach 
molecule experiments, is extremely sensitive to the variation of the electron-to-proton mass 



ratio /i = m e /m p or X e = m e / 'Aqcd [69J- 



S *K S -± = K S * (16) 

A jj, A e 



where K is the enhancement factor. For example, for Cs-Cs collisions we obtained K ~ 400. 
With the Feshbach resonance, however, one is given the flexibility to adjust position of the 
resonance using external fields. Near a narrow magnetic or an optical Feshbach resonance 
the enhancement factor K may be increased by many orders of magnitude. 



IX. MOLECULAR SPECTRA 

Recently we wrote a review about search 
in quasar and laboratory molecular spectra 
to our works. 



or the variation of the fundamental constants 



3- 



Below I present several examples related 



A. Comparison of hydrogen hyperfine and molecular rotational quasar spectra 



The frequency of the hydrogenic hyperfine line is proportional to a 2 fig p atomic units, 

to /i atomic units. The comparison places 
711 ] . Recently similar analysis was repeated 



molecular rotational frequencies are proportiona 
limit on the variation of the parameter F = o?l 

by Murphy et al [72j using more accurate data for the same object at z = 0.247 and for a 
more distant object at z — 0.6847, and the following limits were obtained: 



SF 

— = (-2.0 ±4.4) x 10~ 6 
F y ' 



5F 
~F 



^-1.6 ±5.4) x 10 



-6 



(17) 



(18) 



The object at z = 0.6847 is associated with the gravitational lens toward quasar B0218+357 
and corresponds to the backward time ~ 6.5 Gyr. 
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B. Enhancement of variation of \i in inversion spectrum of ammonia and limit 
from quasar spectra 



Few years ago van Veldhoven et al suggested to use decelerated molecular beam of ND 3 
to search for the variation of jj, in laboratory experiments 73[. Ammonia molecule has a 
pyramidal shape and the inversion frequency depends on the exponentially small tunneling 
of three hydrogens (or deuteriums) through the potential barrier. Because of that, it is 
very sensitive to any changes of the parameters of the system, particularly to the reduced 



mass for this vibrational mode. This fact was used in 28J| to place the best limit on the 
cosmological variation of \x. 

The inversion vibrational mode of ammonia is described by a double well potential with 
first two vibrational levels lying below the barrier. Because of the tunneling, these two levels 
are split in inversion doublets. The lower doublet corresponds to the wavelength A ~ 1.25 cm 
and is used in ammonia masers. Molecular rotation leads to the centrifugal distortion of the 
potential curve. Because of that, the inversion splitting depends on the rotational angular 
momentum J and its projection on the molecular symmetry axis K: 

u inv (J, K) = < v - Cl [J(J + 1) - K 2 } + c 2 K 2 + ■■■, (19) 

where we omitted terms with higher powers of J and K. Numerically, u; mv ~ 23.787 GHz, 
ci » 151.3 MHz, and c 2 « 59.7 MHz. 

In addition to the rotational structure (fT9]) the inversion spectrum includes much smaller 
hyperfine structure. For the main nitrogen isotope 14 N, the hyperfine structure is dominated 
by the electric quadrupole interaction (~ 1 MHz). Because of the dipole selection rule 
AK = the levels with J = K are metastable. In astrophysics the lines with J = K are 
also narrower and stronger than others, but the hyperfine structure for spectra with high 
redshifts is still not resolved. We obtained the following results for NH3 [28j (in atomic 
units): 

^-4.46^. (20) 

^inv ^ 

^ = 5.1^. (21) 

C l,2 A* 

For ND3 the inversion frequency is 15 times smaller and this leads to a higher realative 
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sensitivity of the inversion frequency to fi: 

-cT^s.rf. (22) 

^inv f 1 

-^± = 6.2—. 23 

Cl,2 \i 

We see that the inversion frequency uj® nv and the rotational intervals u inv (Ji, Ki) — 
cJi nv (J2, K2) have different dependencies on the constant \i. In principle, this allows one 
to study time-variation of fi by comparing different intervals in the inversion spectrum of 
ammonia. For example, if we compare the rotational interval to the inversion frequency, 
then Eqs. flU and (J2U) give: 

<Kknv(./l, K l) - ^invih, ^2)]M° nv } = Q g 5/£ 

The relative effects are substantially larger if we compare the inversion transitions with the 
transitions between the quadrupole and magnetic hyperfme components. However, in prac- 
tice this method will not work because of the smallness of the hyperfme structure compared 
to typical line widths in astrophysics. 

We compared the inversion spectrum of NH3 with rotational spectra of other molecules, 
where 

5u^ = 8 ± (25) 

High precision data on the redshifts of NH 3 inversion lines exist for already mentioned object 
B0218+357 at z ~ 0.6847 Comparing them with the redshifts of rotational lines of CO, 
HCO + , and HCN molecules from Ref. |75| one can get the following limit: 

^ = ^ = (-0.6±1.9)xl(r 6 . (26) 

Taking into account that the redshift z ~ 0.68 for the object B0218+357 corresponds to the 
backward time about 6.5 Gyr and assuming linear time dependence , this limit translates 



into the most stringent present limit for the variation rate fi/fi and X e 28]: 



ii/li = X e /X e = (1 ± 3) x 10- 16 yr- 1 . (27) 



A combination of this result and the atomic clock results [2i 
variation of a: 



30| gives the best limit on 



a/a = (-0.8 ±0.8) x 10~ 16 yr" 1 . (28) 
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X. PROPOSALS OF ENHANCED EFFECTS IN DIATOMIC MOLECULES 



In transitions between very close narrow levels of different nature in diatomic molecules 
the relative effects of the variation may be enhanced by several orders of magnitude. Such 
levels may occur due to cancelation between either hyperfine and rotational structures 0] , 
or between the fine and vibrational structures of the electronic ground state 77|. The 
intervals between the levels are conveniently located in microwave frequency range and the 
level widths are very small, typically ~ 10~ 2 Hz. 



A. Molecules with cancelation between hyperfine structure and rotational inter- 
vals 



Consider diatomic molecules with unpaired electron and ground state 2 £. It can be, for 
example, LaS, LaO, LuS, LuO, YbF, etc. 78]. Hyperfine interval Ahf s is proportional to 
a 2 ZF re i(aZ)fig nnc , where F re \ is additional relativistic (Casimir) factor. Rotational interval 
Arot ~ fj, is roughly independent on a. If we find molecule with A^fs ~ A rot the splitting uj 
between hyperfine and rotational levels will depend on the following combination 



uj ~ 



a 2 F Te i(aZ) g nnc - const 



(29) 



Relative variation is then given by 



8uj A 



hfs 



UJ 



UJ 



{2 + K) ^ + ^ 
Q fi'nuc 



(30) 



where factor K comes from variation of F xe \{aZ), and for Z ~ 50, K w 1. Using data from 
Q one can find that uj = (0.002±0.01) cm' 1 for 139 La 32 S 0. Note that for uj = 0.002 cm" 1 
the relative frequency shift is: 



5oj 



UJ 



600 



5a 



a 



(31) 



B. Molecules with cancelation between fine structure and vibrational intervals 



The fine structure interval ujf rapidly grows with nuclear charge Z: 



ry2 2 

ujf ~ Z a , 



(32) 
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The vibration energy quantum decreases with the atomic mass: 

uj vih ~ M" 1 V /2 , (33) 

where the reduced mass for the molecular vibration is M r m p . Therefore, we obtain equation 
Z = Z(M r , v) for the lines on the plane Z, M r , where we can expect approximate cancelation 
between the fine structure and vibrational intervals: 

to — ujf — uovib ~ , v — 1,2, ... (34) 

Using Eqs. fl3"2H3"4"|) it is easy to find dependence of the transition frequency on the funda- 
mental constants: 

5uj if 5a v 5u\ rr / 8a I5a\ , . 

uj uj \ a 2 ii J \ a 2 fi J 

where the enhancement factor K — — determines the relative frequency shift for the given 
change of fundamental constants. Large values of factor K hint at potentially favorable 
cases for making experiment, because it is usually preferable to have larger relative shifts. 
However, there is no strict rule that larger K is always better. In some cases, such as 
very close levels, this factor may become irrelevant. Thus, it is also important to consider 
the absolute values of the shifts and compare them to the linewidths of the corresponding 
transitions. 

Because the number of molecules is finite we can not have u = exactly. However, a 
large number of molecules have uj/ujf 1 and \K\ 3> 1. Moreover, an additional "fine 
tuning" may be achieved by selection of isotopes and rotational, fi-doublet, and hyperfine 
components. Therefore, we have two large manifolds, the first one is build on the electron 
fine structure excited state and the second one is build on the vibrational excited state. If 
these manifolds overlap one may select two or more transitions with different signs of uj. In 
this case expected sign of the variation must be different (since the variation 8oj has the 
same sign) and one can eliminate some systematic effects. Such control of systematic effects 



was used in 



59 



601 ] for transitions between close levels in two dysprosium isotopes. The 
sign of energy difference between two levels belonging to different electron configurations 
was different in 163 Dy and 162 Dy. 

Among the interesting molecules where the ground state is split in two fine structure 
levels and fl34l) is approximately fulfilled, there are Cl^ (enhancement K = 1600), SiBr 
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(K = 360), CuS (K = 24) and IrC (K = 160). The list of molecules is not complete because 



of the lack of data in 78] . The molecules ClJ and SiBr are particularly interesting. For 
both of them the frequency u defined by (j34p is of the order of 1 cm -1 and comparable to 
the rotational constant B. That means that uj can be reduced further by the proper choice 
of isotopes, rotational quantum number J and hyperfine components, so we can expect 
K ~ 10 3 — 10 5 . New dedicated measurements are needed to determined exact values of the 
transition frequencies and find the best transitions. However, it is easy to find necessary 
accuracy of the frequency shift measurements. According to (135]) the expected frequency 
shift is 

= + (36) 

Assuming 5a /a ~ 10~ 15 and uoj ~ 500 cm -1 , we obtain 5uo ~ 10~ 12 cm -1 ~3x 10~ 2 Hz (in 
order to obtain similar sensitivity comparing hyperfine transition frequencies for Cs and Rb 
one has to measure the shift ~ 10~ 5 Hz). This shift is larger than the natural width ~ 10 -2 
Hz. 

C. Molecular ion HfF+ 

The ion HfF + and other similar ions are considered by Cornell's group in JILA for the 
experiment to search for the electric dipole moment (EDM) of the electron. Recent calcula- 
tion by {73] suggests that the ground state of this ion is 1 S + and the first excited state 3 Ai 
lies only 1633 cm -1 higher. Calculated vibrational frequencies for these two states are 790 
and 746 cm -1 respectively. For these parameters the vibrational level v = 3 of the ground 
state is only 10 cm -1 apart from the v — 1 level of the state 3 Ai. Thus, instead of ( |34l) we 
now have: 

" = <** + - '/I « , (37) 
where superscripts and 1 correspond to the ground and excited electronic states. Electronic 
transition uo c \ is not a fine structure transition and (|32l) is not applicable. Instead, we can 
write: 

tUei = uj e \ t o + qx , x = c?ja\ — \. (38) 



Our estimate is [77j] 



5uj ( 2q 5a u e i5u\ ( 5a 5 11 
— w — — + — — w 2000— + 80— 
uj \oj a 2u n J \ a \i / 



(39) 
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5uo « 20000 cm" 1 {8a/ a + 0.045/x///) . (40) 
Assuming ~ 10" 15 we obtain <5c<j ~ 0.6 Hz. The natural width is about 2 Hz. 



We also present the resu 



in 



80 



811 ] . Our estimate is 



t for transition between close levels in Cs2 molecule suggested 



5oo « (-240— - leOO^cm" 1 , (41) 
a ji 



XI. CHANGING PHYSICS NEAR MASSIVE BODIES 



In this section I follow Ref. 82| 



The reason gravity is so important at large scales is that its effect is additive. The same 
should be true for massless (or very light) scalars: its effect near large body is proportional 
to the number of particles in it. 

For not-too-relativistic objects, like the usual stars or planets, both their total mass M 
and the total scalar charge Q are simply proportional to the number of nucleons in them, 
and thus the scalar field is simply proportional to the gravitational potential 

cf) - O = K(GM/rc 2 ) . (42) 

Therefore, we expect that the fundamental constants would also depend on the position via 
the gravitational potential at the the measurement point. 

Gravitational potential on Earth is changing due to ellipticity of its orbit: the correspond- 
ing variation of the Sun graviational potential is 5(GM/rc 2 ) = 3.3 • 10~ 10 . The accuracy of 
atomic clocks in laboratory conditions is about 10 -16 . As an example we consider recent 
work 30j who obtained the following value for the half-year variation of the frequency ra- 
tio of two atomic clocks: (i) optical transitions in mercury ions 199 Hg + and (ii) hyperfine 
splitting in 133 Cs (the frequency standard). The limit obtained is 

«n(^) = (0.7±1.2)-10- 15 (43) 

^Cs 

For Cs/Hg frequency ratio of these clocks the dependence on the fundamental constants was 
evaluated in [29(] with the result 

6ln e*i) = - o,o /K/ A ^) _ *K/Mp) (44) 

u Cs « {m q /A QC D) {m e /M p ) 
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Another work 83j compare H and 133 Cs hyperfine transitions. The amplitude of the half- 
year variation found were 

\5ln{u H /u Cs )\<l-lQ- 15 (45) 



The sensitivity 29j 



8ln{ 



-0.83^-0.11^/^ 



(46) 



'u Cs ' a (m q /A QC D) 

There is no sensitivity to m e /M p because they are both hyperfine transitions. 

As motivated above, we assume that scalar and gravitational potentials are proportional 

to each other, and thus introduce parameters fcj as follows 



5a , r ,GM, 



a 



rcr 



S(m q /A QCD ) ( GM 
(m q /A QOD ) q [ re- ] 

6(m e /A QCD ) = 6(m e /M p ) = fc ^GM, 
{m e /A QCD ) (m e /M p ) e rc 2 ' 



(47) 

(48) 
(49) 



where in the r.h.s. stands half-year variation of Sun's gravitational potential on Earth. 
In such terms, the results of Cs/Hg frequency ratio measurement 30] can be rewritten as 



k a + 0.17k e = (-3.5 ±6) ■ 10 



-7 



(50) 



The results of Cs/H frequency ratio measurement [83j can be presented as 



\k a + 0.13*, | < 2.5 • 10" 



(51) 



Finally, the result of recent measurement 84| of Cs/H frequency ratio can be presented as 



k a + 0.13fc g = (-1 ± 17) • 10 



-7 



(52) 



The sensitivity coefficients for other clocks have been discussed above. 

Two new results have been obtained recently. From transition between close levels in Dy 
we obtained js^] 

k a = (-8.7 ±6.6) ■ 10~ 6 (53) 



From optical Sr/hyperfine Cs comparison we obtained 86] 



k a + 0.36k e = (1.8 ±3.2) ■ 10" 



(54) 
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Combination of the data gives 86] 



k Q = (-2.3 ±3.1) ■ l(T fc 

k e = (1.1 ± 1.7) • 10 -5 
k a = (1.7 ±2.7) • 1(T 5 



(55) 

(56) 
(57) 
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